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The standard engineering approach to modelling of complex systems is 
highly compositional. In order to be able to understand (or to control) the 
behavior of a complex dynamical systems, it is often desirable, if not necessary, 
to view this system as an interconnection of smaller interacting subsystems, 
each of these subsystems having its own functionalities. In this paper, we pro- 
pose a compositional approach to the computation of bisimulation functions 
for dynamical systems. Bisimulation functions are quantitative generaliza- 
tions of the classical bisimulation relations. They have been shown useful for 
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I/*-) ■ simulation-based verification or for the computation of approximate symbolic 



abstractions of dynamical systems. In this technical note, we present a con- 
structive result for the composition of bisimulation functions. For a complex 



dynamical system consisting of several interconnected subsystems, it allows 
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us to compute a bisimulation function from the knowledge of a bisimulation 



function for each of the subsystem. 
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1 Introduction 



The standard engineering approach to modelling of complex systems is highly com- 
positional. In order to be able to understand (or to control) the behavior of a com- 
plex dynamical systems, it is often desirable, if not necessary, to view this system 
as an interconnection of smaller interacting subsystems, each of these subsystems 
having its own functionalities. System on chips, for instance, are often complex 
circuits that can be decomposed into smaller (and thus simpler) circuits. 

Albeit the simplification of the modelling process, a modular representation of 
complex systems can greatly simplify the analysis process. In computer science, 
compositionallity and concurrency [13] have been a very active research field. In 
the system engineering science, a compositional approach is also often used (see 
e.g. [10J). In this paper, we propose a compositional approach to the computation 
of bisimulation functions for dynamical systems. 

Bisimulation functions have been introduced in [5] as a quantitative general- 
ization of the classical notion of bisimulation relations that have been extensively 
and successfully used in purely discrete systems analysis [2] . Bisimulation functions 
measure how far two states of a system are from being bisimilar, thus enabling the 
quantification of the distance between trajectories originating from different states. 
Thus, these functions allow us to define a natural notion of neighborhood for trajec- 
tories of a system. Recently, several promising papers have shown that bisimulation 
functions can be used for simulation-based verification [SJ HJ (III H2] or for the com- 
putation of approximate symbolic abstractions of dynamical systems [TJ [H] . 

In this technical note, we present a constructive result for the composition of 
bisimulation functions. For a complex dynamical system consisting of several in- 
terconnected subsystems, it allows us to compute a bisimulation function from the 
knowledge of a bisimulation function for each of the subsystem. Similar to Lya- 
punov functions for interconnected systems [10], a small gain condition has to be 
fulfilled in order be able to compose bisimulation functions. The paper is organized 
as follows. First, we present the notion of interconnection of subsystems useful for 
compositional modelling of dynamical systems. Then, we introduce the notion of 
bisimulation function and develop a result on composition of bisimulation functions. 



2 Compositional Modelling of Dynamical Systems 

Compositional modelling allows us to see a complex dynamical system S as a set of 
several smaller subsystems Si, ... , S m , interacting together. This is a standard en- 
gineering approach and softwares such as Simulink or Scicos gained their popularity 
from the possibility of modular representation of complex systems. In this section, 
we present the notion of interconnection of subsystems useful for compositional 
modelling of dynamical systems. In the following, we only define the interconnec- 
tion of two subsystems; however, the extension to systems with more components is 
straightforward (see e.g. [TB]). 

Let us consider two dynamical systems, Si and S2 of the following form: 

Si : &i(t) = fi(xi(t),Ui(t)), i = 1,2. 

where X{(t) G M. ni and ««(£) G M mi denote the state and input variables of Sj. The 
input vector is of the form Ui(t) = [vi(t),Wi(t)], where t>j(£) G M Pi denotes the inputs 
used for the interconnection of Si and S 2 and Wi(t) G M. qi denotes the external 
inputs (see Figured]). 
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Figure 1: Subsystem Sj. 



The interconnection of S x and S 2 is achieved by feeding the system inputs Vi(t) 
and v 2 (t) with the state variables x 2 (t) and Xi(t) (see Figure [2]). We therefore assume 
that pi = n 2 and p 2 = n\. Then, the interconnection of Si and S 2 is formally defined 
as follows: 

Definition 1 The interconnection of Si and S 2 is the dynamical system S given 
by the differential equation 

= fi(xi(t),x 2 (t),Wx(t)), 
= f2(x 2 (t),x 1 (t),w 2 (t)) 
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Figure 2: The composition of Si and S 2 . 

The state of S is x(t) = [xi(t), X2{t)} G M™ with n = n\ + n 2 and the input of S is 
u(t) = [wi(t), w 2 (t)} G M m with m = qi + q 2 . Then, the system E can be written 
under the form 

Z:x(t) = f(x(t),u(t)), 

that is similar to Si and S 2 . Then, this means that S can be composed with another 
system, enabling the hierarchical modelling of dynamical systems. 

3 Composition of Bisimulation Functions 

We first present the notion of bisimulation function, then we will give a result on 
composition of bisimulation functions. 



3.1 Bisimulation functions 

Let us consider a dynamical system of the form 

Z:x{t) = f{x{t),u{t)) 

where x(t) G MJ 1 and u(t) G M. m . Bisimulation functions have been introduced in [5] 
as a quantitative generalization of the classical notion of bisimulation relations that 
have been extensively and successfully used in purely discrete systems analysis [2]. 



Bisimulation functions measure how far two states of a system are from being bisim- 
ilar, thus enabling the quantification of the distance between trajectories originating 
from different states. Thus, these functions allow us to define a natural notion of 
neighborhood for trajectories of a system. The following definition slightly differs 
from the original definition in [5] . It is the continuous time version of the definition 
given in [I] which makes it suitable for simulation-based verification. 

Definition 2 A smooth function V : M. n xR"4 R + is a bisimulation function for 
Eif 

\\x — x'\\ < V(x, x) (1) 

and there exists A > 0, 7 > such that for all x E R n , x' E R n , Vu E JR m , u' E W n , 

dV dV 

—f(x, u) + —f(x', u') < -XV(x, x') + 1 \\u- u'\\ . (2) 

Bisimulation functions have the following property which makes them suitable 
tools for simulation-based verification [HI IH dH [12] or for the computation of ap- 
proximate symbolic abstractions of dynamical systems [7J Q3] . 

Theorem 1 Let us consider x{t) and x'{t) be the trajectories ofT, given by 

x(t) = f(x(t),u(t)) and x' (t) = f(x'{t),u\t)). 

Then, we have for all t > 

\\x(t) - x'(t)\\ < V(x(t), x'(t)) < e- xt V(x(0),x'(0)) + ^\\u - u'W^ 

A 

where \\u — u'\\oo = sup 4>0 \\u(t) — u'(t)\\. 

Proof: From equation ([T]), we have the first inequality. From equation (T5]), we 
have 

dV(x(t),x'(t)) 



dt 



< -AV r (x(t),x / (t))+7||w(t)-M , (t)|| < -XV{x{t),x'{t)) + -f\\u-u'\ 



Let r](t) = e ^^(^(O), x'(0)) + ?||w — u'\\oo, it is a solution of the differential equation 

rj(t) = -\rj(t) + 7||w - w'Hoo. 

Moreover, V(x(0), x'(0)) < n(0); then, from the funnel theorem (see e.g. (9]), it 
follows that for all t > 0, V{x{t), x'{t)) < r)(t). ■ 



The practical computation of bisimulation functions is out of the scope of this 
technical note. However, we refer the interested reader to [21 E] for computational 
methods applying to linear and nonlinear dynamical systems. 

Let us remark that the previous theorem clearly shows the existing relation 
between the notion of bisimulation function and the notion of incremental input-to- 
state stability pQ (close initial states and close inputs lead to close trajectories of E). 
This connection was already pointed out in the work [14] where incremental input- 
to-state stability was shown sufficient for the existence of approximately bisimilar 
symbolic abstractions of a dynamical system. 

3.2 A Composition Result for Bisimulation Functions 

We now consider the problem of composing bisimulation functions. For complex sys- 
tems that consists of several interconnected subsystems, it is interesting to develop 
compositional analysis methods. Let us assume that we are given a bisimulation 
function for each subsystem, then the question is whether it is possible or not to 
compose these functions to design a bisimulation function for the global system. The 
following result shows that the composition is possible under a small gain condition. 
It has similarities with [10] where a compositional result for ISS-Lyapunov functions 
is developped. 

Theorem 2 Let Ei and E2 be dynamical systems and let E be the interconnection 
0/E1 and E 2 as defined in Definition 1. Let V± and V2 be simulation functions for Ei 
and E 2; we denote by Ai and j% (respectively A 2 and 72 ) the real numbers such that 
equation (TJ|) holds for V\ (respectively V2). Then, under the small gain condition 
YY~ < 1, there exists V a bisimulation function for E of the form: 

V(x,x') = aiVi(xi,z' 1 ) + ct2V2(x2,x 2 ) where x = [rci,^], x' = [x^^x^]. (3) 

The couple (ai,a 2 ) can be chosen as follows 

g < «i < ^ and a 2 = l if A x < 72 

«i = l and f 2 <a 2 <± tfA 2 <7i (4) 

ax — 1 and «2 = 1 in the other cases. 

Proof. Let V be a function of the form ([3]), we look for conditions on a\ and ct2 
ensuring that V is a bisimulation function for E. First, let us remark that if a± > 1 
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and a 2 > 1 then, 

V(x, x') > Vi(Xl, x'i) + V2(X2,x' 2 ) > \\Xi — x\\\ + ||a;2 — X 2 \ 

because V\ and V2 satisfy equation flTJ. Then, by remarking that 



|| X — X || — v/ \\X\ — ^l|| ~r r-^2 — *^2 II — 11^1 — *^1 + 11*^2 — ^9 ' 

it follows that V satisfies equation ([1]) as well. Let u = [u>i,u> 2 ], u' = [w[,w 2 ] be 
inputs of E. Then, we have 

— f{x,U) + —f{x,U) = a 1 —f 1 {x 1 ,X2,W 1 )+a 2 -K-f2{X2,Xl,W2) 

+01-0-7/1 (»1, X 2 , Wi) + OL 2 — j 2 [X 2 , X 1 ,W 2 ) 

< «i {-X\V 1 {xx,x' 1 ) + 7i||[a;2,wi] - [x 2 ! w i]||) 

+a 2 (-A 2 V2(^2,4) + «272 1| [2:1,^2] - [a^mil) 

because V\ and V2 satisfy equation (J2J). Further, we have 

||[x 2 ,u'i] - [x 2 ,wi]|| = a/||x 2 - x' 2 || 2 + II wi - w[\\ 2 < \\x 2 - x 2 \\ + ||wi - m| 
and 



|| [x\, w 2 ] — [x[, w 2 ]\\ = vjjxi — x'JI 2 + ||w 2 — if 2 || 2 < ||xi — aJiH + ||u> 2 — w' 2 \\. 

Therefore, 

dV dV 

—f(x,u) + —f(x',u) < ai(-XiVi(xx,x x ) +7i||x 2 - z 2 ll +7i||wi -will) 

+a 2 (-A 2 V 2 (a; 2 ,x 2 ) + 7 2 ||xi - x[\\ + 7 2 ||w 2 - w' 2 \\ 
Then, since V\ and V 2 satisfy equation (jTJ), it follows that 

—f(x,u) + —f(x',u) < a.\ {-\\V\{x\,x' x ) + 7iV 2 (x 2 ,x 2 ) +7i||wi - w[\\) 

+a 2 (-A 2 V2(x 2 ,a; 2 ) + 7 2 Vi(a;i,x / 1 ) + 7 2 ||w 2 - w 2 | 
< -(ckiAi - a 2 ^ 2 )Vi(xi,x[) H-ai7i||wi - w[\\ 
-(a 2 A 2 -a 1 ji)V 2 (x2,x' 2 ) + a 2 7 2 ||u> 2 - w' 2 \\. 



Let us assume that ctiAi — a 2 7 2 > and a 2 A 2 — «i7i > 0, then let us define 

. .'aiAi -a 2 72 a 2 A 2 -ai7i\ 
A = mm I , I and 7 = «i7i + a 2 7 2 . 



«i a 2 

By remarking that \\wi — w[\\ < \\u — u'\\ and ||w 2 — w' 2 \\ < ||« — tt'|| it follows that 

9y dV 

— f(x, u) + Q-;f(z, u) < -XV(x, x) + j\\u - u'|| . 

Therefore, we proved that if a\ > 1, a 2 > 1, «iAi — a 2 7 2 > and a 2 A 2 — a\ji > 0, 
then V is a bisimulation function for S. Let us show that a necessary and sufficient 
condition for the existence of a couple (ai,a 2 ) satisfying these four inequalities is 
YY- < 1. Let the inequalities hold, then particularly, 

ctiAiAi > a 2 7 2 A 2 > aijij 2 - 



It follows that necessarily 7^~ < 1. Conversely, if ^y- < 1, there are only three 
possible configurations shown on Figures El @] and Then, by choosing a x and a 2 
as in equation (jlj), the four inequalities hold. ■ 



a 2 ii 




«i 



Figure 3: Configuration 1: Xi < 72 . 

This theorem provides us with a method to compute compositionally bisimu- 
lation functions for composite systems. Note that it is subject to a small gain 
condition that is y~\~ < 1- Let us remark that the choice of the couple (ai,a 2 ) 
given in equation (jlj) is only one possible choice among many others satisfying the 
inequalities «i > 1, a 2 > 1, ctiAi — a 2 7 2 > and a 2 A 2 — «i7i > 0. Another choice 
might be more suitable depending on the considered application. 
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72«2 = Ai«i 
«2 A 




\ 2 a 2 = 71CK1 



«2 = 1 



Q-4 



Figure 4: Configuration 2: A2 < 71 ■ 
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Figure 5: Configuration 3: other cases. 

4 Conclusions 

In this technical note, we presented a constructive result for the composition of 
bisimulation functions. For a complex dynamical system consisting of several in- 
terconnected subsystems, it allows us to compute a bisimulation function from the 
knowledge of a bisimulation function for each of the subsystem. Similar to Lyapunov 
functions for interconnected systems, a small gain condition has to be fulfilled in 
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order be able to compose bisimulation functions. 

In the context of the VAL-AMS project, this result shall be useful for the com- 
putation of bisimulation functions for large scale analog circuits which can be seen 
as the interconnection of smaller circuits. The knowledge of a bisimulation function 
is required for simulation-based verification [51 HI [UJ [12] or for the computation of 
approximate symbolic abstractions of dynamical systems [TJ IT4"] . 

References 

[1] D. Angeli. A Lyapunov approach to incremental stability properties. IEEE 
Trans. Automatic Control, 47(3):410-421, March 2002. 

[2] E.M. Clarke, O. Grumberg, and D. Peled. Model Checking. MIT Press, 1999. 

[3] A. Girard. Reachability of uncertian linear systems using zonotopes. In Hybrid 
Systems : Computation and Control, volume 3414 of LNCS, pages 291-305. 
Springer, 2005. 

[4] A. Girard. Simulation-based techniques for verification of dynamical systems. 
Technical report, 2007. VAL-AMS Deliverable. 

[5] A. Girard and G. J. Pappas. Approximation metrics for discrete and continuous 
systems. IEEE Transactions on Automatic Control, 52 (5): 782-798, May 2007. 

[6] A. Girard and G.J. Pappas. Approximate bisimulation relations for constrained 
linear systems. Automatica, 43(8), 2007. 

[7] Antoine Girard. Approximately bisimilar finite abstractions of stable linear 
systems. In Hybrid Systems: Computation and Control, volume 4416 of LNCS, 
pages 231-244. Springer, 2007. 

[8] Antoine Girard and George J. Pappas. Verification using simulation. In Hybrid 
Systems: Computation and Control, volume 3927 of LNCS, pages 272-286. 
Springer, 2006. 

[9] J.H. Hubbard and B.H. West. Differential equations: a dynamical systems 
approach. Springer, 1995. 

10 



[10] Z.P. Jiang, I.M.Y. Mareels, and Y. Wang. A lyapunov formulation of nonlinear 
small gain theorem for interconnected iss systems. Automatica, 32(9):1211- 
1215, 1996. 

[11] A. Agung Julius, Georgios E. Fainekos, Madhukar Anand, Insup Lee, and 
George J. Pappas. Robust test generation and coverage for hybrid systems. 
In Hybrid Systems: Computation and Control, volume 4416 of LNCS, pages 
329-342. Springer, 2007. 

[12] F. Lerda, J. Kapinski, E.D. Clarke, and B.H. Krogh. Verification of supervi- 
sory control software using state proximity and merging. In Hybrid Systems : 
Computation and Control, LNCS. Springer, 2008. To appear. 

[13] R. Milner. Communication and concurrency. Prentice Hall, 1989. 

[14] G. Pola, A. Girard, and P. Tabuada. Approximately bisimilar symbolic models 
for nonlinear control systems. 2007. Submitted. 

[15] Y. Tazaki and J. Imura. Bisimilar finite abstractions of interconnected systems. 
In Hybrid Systems : Computation and Control, LNCS. Springer, 2008. To 
appear. 



11 



